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Abstract. A derivation operator and a divergence operator are defined on 
the algebra of bounded operators on the symmetric Fock space over the com- 
plexification of a real Hilbert space t) and it is shown that they satisfy simi- 
lar properties as the derivation and divergence operator on the Wiener space 
over f). The derivation operator is then used to give sufficient conditions for 
the existence of smooth Wigner densities for pairs of operators satisfying the 
canonical commutation relations. For tj = the divergence operator is 

shown to coincide with the Hudson-Parthasarathy quantum stochastic integral 
for adapted integrable processes and with the non-causal quantum stochastic 
integrals defined by Lindsay and Belavkin for integrable processes. 



1. Introduction 

Infinite-dimensional analysis has a long history: it began in the sixties (work 
of Gross ||Gro67|| , Hida, Elworthy, Kree, ...), but it is Malhavin ||Mal78|| who 



has applied it to diffusions in order to give a probabilistic proof of Hormander's 
theorem. Malliavin's approach needs a heavy functional analysis apparatus, as 
the Ornstein-Uhlenbeck operator and the definition of suitable Sobolev spaces, 
where the diffusions belong. Bismut ||Bis81|| has given a simpler approach based 



upon a suitable choice of the Girsanov formula, which gives quasi-invariance 
formulas. These are differentiated, in order to get integration by parts formulas 
for the diffusions, which where got by Malliavin in another way. 

Our goal is to generalize the hypoellipticity result of Malliavin for non-com- 



mutative quantum processes, by using Bismut's method, see also |[FLS99|| . For 
that we consider the case of a non-commutative Gaussian process, which is the 
couple of the position and momentum Brownian motions on Fock space, and we 
consider the vacuum state. We get an algebraic Girsanov formula, which allows 
to get integration by parts formulas for the Wigner densities associated to the 
non-commutative processes, when we differentiate. This allows us to show that 
the Wigner functional has a density which belongs to all Sobolev spaces over M^. 
Let us remark that in general the density is not positive. 

If we consider the deterministic elements of the underlying Hilbert space of the 
Fock space, the derivation of the Girsanov formula leads to a gradient operator 
satisfying some integration by parts formulas. This shows it is closable as it 
is in classical infinite-dimensional analysis. But in classical infinite-dimensional 
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analysis, especially in order to study the Malliavin matrix of a functional, we need 
to be able to take the derivation along a random element of the Cameron-Martin 
space. In the commutative set-up, this does not pose any problem. Here, we 
have some difficulty, which leads to the definition of a right-sided and a left-sided 
gradient, which can be combined to a two-sided gradient. 

We can define a divergence operator as a kind of adjoint of the two-sided 
gradient for cylindrical (non-commutative) vector fields, but since the vacuum 
state does not define a Hilbert space, it is more difficult to extend it to general 
(non-commutative) vector fields. 

We show that the non-commutative differential calculus contains in some sense 
the commutative differential calculus. 

In the white noise case, i.e. if the underlying Hilbert space is the L^-space of 
some measure space, the classical divergence operator defines an anticipating sto- 
chastic integral, known as the Hitsuda-Skorohod integral. We compute the matrix 
elements between exponential vectors for our divergence operator and use them 
to show that the divergence operator coincides with the non-causal creation and 
annihilation integrals defined by Belavkin |[Bel91a|, [Bel91b|| and Lindsay |[Lin93 



for integrable processes, and therefore with the Hudson-Parthasarathy [pP84 
integral for adapted processes. 

2. Analysis on Wiener space 

Let us first briefiy recall a few definitions and facts from analysis on Wiener 
space, for more details see, e.g., ||Jan97| , |Mal97| , [Nua95| , [Nua98| , |Ust95|| . Let f) be a 



real separable Hilbert space. Then there exists a probability space (fi, JF, P) and 
a linear map W : i) ^ L^(f2) such that the W{h) are centered Gaussian random 
variables with covariances given by 

E{W{h)W{k)) = {h,k), for all /i,A; G f^. 

Set Hi = W{i)), this is a closed Gaussian subspace of L'^{Q) and 1^ : f) — >• Tii C 
L'^{Q) is an isometry. We will assume that the a-algebra JF is generated by the 
elements of Tii. We introduce the algebra of bounded smooth functionals 

S = {F = f{W{h), W{hr,)) \neN,fe cr(M"), h,...,Ke f)}, 

and define the derivation operator D : S L'^{Q) ® f) = L^(i7; [)) by 

i=l * 

for F = f(W{hi), . . . , W{hn)^ G S. Then one can verify the following properties 
of £). 

1. D is a derivation (w.r.t. the natural L°°(n)-bimodule structure of L'^{Q] [))), 
i.e. 

D{FG) = F{DG) + {DG)F, for all F,GeS. 
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2. The scalar product {h, DF) coincides with the Prechet derivative 

d 



/(W-'Ci) + Ai), . . . , W{K) + E{h, K)) 

£=0 

for all F = f{W{hi), ... , W{hn)) e S and all h e i). 

3. We have the following integration by parts formulas, 

(2.1) E{FW{h)) = E{{h,DF)) 

(2.2) E{FGW{h)) = E{{h,DF)G + F{h,DG)) 

ioT all F,G eS, hei). 

4. The derivation operator D is a closable operator from U'{Q) to L*'(Q; f)) for 
1 < p < oo. We will denote its closure again by D. 

We can also define the gradient D^F = {u, DF) w.r.t. f)- valued random variables 
u G L^{fl; [)), this is L°°(0)-linear in the first argument and a derivation in the 
second, i.e. 

DpuG = FDuG, 
D^{FG) = F{D^G) + {D^F)G. 

L^(f2) and L'^{Q; [)) are Hilbert spaces (with the obvious inner products), 
therefore the closability of D implies that it has an adjoint. We will call the 
adjoint of D : L^{Q) — >• L^{Q; ()) the divergence operator and denote it by 
5 : L\Q; t)) L'^{Q). Denote by 

5(, = jxi = ^ Fj <^hj neN,Fi,... ,FneS,hi,... ,hnei) 
the smooth elementary f)-valued random variables, then S{u) is given by 

n n 

d{u) = Y,F,Wihj) - Y.{hj,DF,) 

for u = E^=i^i ® hj e St,. If we take, e.g., i) = L^{R+), then Bt = M/(l[o,t]) 
is a standard Brownian motion, and the F)-valued random variables can also be 
interpreted as stochastic processes indexed by M_|_. It can be shown that S{u) 
coincides with the Ito integral J^^ UtdWt for adapted intcgrablc processes. In 
this case the divergence operator is also called the Hitsuda- Skorohod integral. 

The derivation operator and the divergence operator satisfy the following re- 
lations 

(2.3) Dhi^i^)) = {h,u) + ~6{Di,u), 

(2.4) E{6{u)6{v)) = E{{u,v)) +E{Tt{Duo Dv)), 

(2.5) S{Fu) = FS{u) - {u,DF), 
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for /i G [), M,f G F G S. Here D is extended in the obvious way to f)-valued 
random variables, i.e. as D^idfj. Thus Du is an {)®f)-valued random variable and 
can also be interpreted as a random variable whose values are (Hilbert-Schmidt) 
operators on f). If {ej;j G N} is a complete orthonormal system on 1), then 
Tr{Du o Dv) can be computed as Tt{Du o Dv) = X]r?=i -^e, (m, ^j)De^ {v, Cj). 

3. The non-commutative Wiener space 

Let again f) be a real separable Hilbert space and let f)c be its complexification. 
Then we can define a conjugation ~ : t)c ^ f)c by /ii + ih2 = hi — ih2 for 
hi,h2 G {)c- This conjugation satisfies {h,k') = {h,k) = {k,h) for all h,k E P)c- 
The elements of f) are characterized by the property h = h, we will call them real. 

Let S) = Ts{[)c) be the symmetric Fock space over i)c, i.e. = 0„gp| t)c"' 
where '©' denotes the symmetric tensor product, and denote the vacuum vector 
1 + + ■ ■ ■ by n. It is well-known that the symmetric Fock space is isomorphic 
to the complexification of the Wiener space L'^{Q) associated to f) in Section ^. 
We will develop a calculus on the non-commutative probability space {B{S}),K), 
where E denotes the state defined by E(X) = {n,Xn) for X G B{S^). To 
emphasize the analogy with the analysis on Wiener space we call {B{S^),K) the 
non-commutative Wiener space over f). 

The exponential vectors {£{k) = Yl'^=o ^ ^ ^'c} total in S), their scalar 
product is given by 

{S{ki),S{k2)) = e^''''^\ 

We can define the operators a{h),a'^{h),Q{h), P{h) (annihilation, creation, 
position, momentum) and U{hi,h2) with h,hi,h2 G P)c on Sj, see, e.g., ||Bia93| , 
Mey95|, [Par92|] . The creation and annihilation operators a~^{h) and a{h) are 
closed, unbounded, mutually adjoint operators. The position and momentum 
operators 

Q{h) = {a(h) + a+(/i)), and P{h) = i{a(h) - a+(/i)) 

are self-adjoint, if h is real. 

The commutation relations of creation, annihilation, position, and momentum 
are 

[a{h),a+{k)] = {h,k), [a{h),a{k)] = [a+{h),a+{k)] = 0, 

[Qih),Qik)] = [Pih),Pik)]=0, [Pih),Qik)] = 2i{h,k). 

The Weyl operators U{hi,h2) can be defined by U{hi,h2) = exp (iP(/ii) -|- 
iQ{h2)) = exp2(a(/i2 — ihi) + a"'"(/i2 — ihi)), they satisfy 

U{hi, h2)U{ki, k2) = expi((/i2, A^i) - {h, k2))U{hi + /i2, ki + ^2) 

Furthermore we have U {hi, h2)* = U{—hi, — /12) and U {hi, /i2)~^ = U{—hi, — /i2)- 
We see that U{hi,h2) is unitary, if hi and /12 are real. These operators act on 
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the vacuum Q = S{0) as 

U{hi, h2)il = exp I j £ {hi + z/ia) 

and on general exponential vectors £ (/) = Yl'^=o 

U{h,, h2)£ if) = exp (^-(7, h + th) - ^J}llM±ihM^ £{f + h+ ih^) . 

The operators a{h),a^{h),Q{h), P{h) and U{hi,h2) are unbounded, but their 
domains contain the exponential vectors. We will want to compose them with 
bounded operators on i^, to do so we adopt the following convention. Let 



L{£{i)c),^) = [Be Lin(span(£({)c)),-^) BB* G Lin (span(£([)c)), ^) 

s.t. {£if),B£ig)) = {B*£{f),£{g)) for all f,g E f)c}, 

i.e. the space of linear operators that are defined on the exponential vectors and 
that have an "adjoint" that is also defined on the exponential vectors. Obviously 
a{h),a'^{h),Q{h), P{h),U{hi,h2) G C[£{[)c),f)). We will say that an expres- 
sion of the form X]j=i -^jBjYj with Xi, . . . , X„, Yi, . . . ,Yn G C(^£{i)c), S^) and 
Bi, . . . ,Bn G B{S)) defines a bounded operator on i^, if there exists a bounded 
operator M G B{S)) such that 

n 

{£{f),M£{g)) = J2{X;£{f),B,Y,£{g)) 

holds for all f,gE P)c- If it exists, this operator is unique, because the exponential 
vectors are total in Sj. We will then write 

n 

i=i 

4. Weyl calculus 
Definition 4.1. Let h = (h, /la) G f) ® M^. We set 

DomO;, = : ^ C 3M G i3(^),VA;i,A;2 G f)c : (^(fci),M£:(fc2)) = 
(4.1) ^ j {£{ki), U{uhi, vh2)£{k2))J'-^v{u, v)dudv^ 

and for ip G DovnOh we define Oh{'-p) to be the bounded operator M appearing in 
Equation (|4.1| ), it is uniquely determined due to the totality of {£{k) : k G Pic}- 



We take the Fourier transform JF as 



T(p{u^ v) = — / (p{x, y) exp {i{ux + vy))(ix(iy. 
2vr J^2 
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Its inverse is simply 

J-'^^(p{x,y) = — / (p{u,v) exp ( — i{ux + vy))dudv. 

Remark 4.2. If is a Schwartz function on M^, then one can check that Oh{^p) = 
2~: /ir2 -^~^V'(w, f ) exp {iuP{hi) + ivQ{h2))du(iv defines a bounded operator. It 
is known that the map from iS(]R^) to B{SS) defined in this way extends to a 
continuous map from L^(R) to B{SS) for all p G [1, 2], but that for p > 2 there 
exist functions in Lp(R^) for which we can not define a bounded operator in this 
way, see, e.g., | |Won98|| and the references cited therein. But it can be extended 
to exponential functions, since ^J-'~^ expi{xQU + y^v) = 5(xo,j/o) thus 

Oh{expi{xou + yov)) = [/(xo/ii, ^0^2)- 

Lemma 4.3. Let I < p < 2 and /i G f) ® such that {hi, ^12) 7^ 0. Then we 
have L^{R'^) C Dom Oh and there exists a constant Ch,p such that 

\\0h{<^)\ \ < Ch,p\\<^\\p 

for allipe Lp(M2). 



Proof. This follows immediately from ||Won98| , Theorem 11.1], where it is stated 



for the irreducible unitary representation with parameter h = 1 of the Heisenberg- 
Weyl group. □ 

As 'joint density' of the pair (P(/ii), (5(^2)) we will use its Wigner distribution. 

Definition 4.4. Let $ be a state on B(Sj). We will call dVF/i$ the Wigner 
distribution of (P(/ii), Q(/i2)) in the state $, if 

is satisfied for all Schwartz functions if. 

In general, dWh,<i> is not positive, but only a signed measure, since Oh does 
not map positive functions to positive operators. But we can show that it has a 
density. 

Proposition 4.5. Let h = (hi, h2) G f) ® ffi^ such that {hi, /i2) 7^ and let $ be 
a state on B{S)). Then there exists a function Wh,<s> £ n2<p<oo -^^(^^) ■^'"^^ ^^'^^ 
dVI4,* = Wh,<!>dxdy. 



Proof. It is sufficient to observe that Lemma ^]3| implies that the map ip 



$(0/i(<y9)) defines a continuous linear functional on L'p(K'^) for 1 < p < 2. □ 

The following proposition will play the role of the Girsanov transformation in 
classical Malhavin calculus. If we conjugate Oh{^p) with f/(— /C2/2, /ci/2) for k G 
[)(8>M^, then this amounts to a translation of the argument of (p by ((fci, hi), {k2, h2)). 
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Proposition 4.6. Let h,k E f) (8) and ip e DomO^j. Then we have 
U{-k2/2,h/2)OMU{-k2/2,h/2r = On{T^^k,M),{k2M))V) 
where ^^^^y^)ip{x, y) = ip{x + xq, y + yo)- 
Proof. For {u,v) e M?, we have 

U{-k2/2, ki/2) exp {i{uP{hi) + vQ{h2))) U{-k2/2, 
= U{-k2/2, kj2)U{uhi,vh2)U{-k2/2, ki/2y 
= e:K.p -i{u{ki, hi) + v{k2,h2))U{uhi,vh2) 
and therefore 

U{-k2/2, ki/2)On{v)U{-k2/2, A;i/2)* 

T~^if{u, v) exp (^—i(u{ki, hi) + v{k2, h2))) expi(uP{hi) + vQ{h2))dudv 
= / J^"^T((fc^,/j^),(fe2,ft2»v3(M,w)expi(MP(/ii) +^;(5(/^,2))d^^d^; 

= Oh {T{{kuhi),{k2M))V) ■ 

□ 

Prom this formula we can derive a kind of integration by parts formula that 
can be used to get the estimates that show the differentiability of the Wigner 
densities. 

Proposition 4.7. Let h e M."^ , k e i)c <S) C^, and ip such that <y^, |f , || e 

DomO/j. Then [Q{ki) — P{k2),0h{f)] defines a bounded operator on S) and we 

have 

'-[Q(ki) - P(k2),OM] = Oh ({ki,hi)^ + {k2M)^ 

Proof. For real k this is the infinitesimal version of the previous proposition, just 
differentiate 

U{ek2/2,eki/2)On{^)U{ek2/2,eki/2y = Oh{T^e{k,M)Ak2M))^) 
with respect to e and set £ = 0. For complex k it follows by hnearity □ 

Like the integration by parts formula in classical Malliavin calculus, this for- 
mula follows from a Girsanov transformation. Furthermore, it can also be used 
to derive sufficient conditions for the existence of smooth densities. 

Proposition 4.8. Let k e N, /i e f) (8) with {hi, /i2) 7^ 0, and $ a vector state, 
i.e. there exists a unit vector uj G S) such that ^{X) — {uj^Xuj) for all X e B{S)). 
If there exists a k E l)c <S) such that 



uje Pi DomQ{kiy'P{k2)^'n Pi DomQ{kiy'P{k2y 



\K2 
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and 



{hi,ki)^0 and (/i2, ^2) 7^ 0, 
HP''^(M?), i.e. the Wigner density Wh,^ lies in the Soholev 



then Wh,<i, e n2<p<oo 
spaces of order n for all 2 < p < 00 . 

Proof. We will show the result for k, = 1, the general case can be shown similarly 
(see also the proof of Theorem |7l2|) . Let iphe a Schwartz function. Let p G [1, 2]. 
Then we have 



ox 



ox 



< 



Similarly, we get 



dip 
dy 



dW, 



< 



Ch,p{\mkiM + \\Q(kiM) 

2\{h,h,)\ 
C,,J||P(A;2V|| + ||P(fc2VII),, ,, 



\p- 



2\{h,h2) 



\pi 



and together these two inequalities imply Wh,<!> G H^' 



,l/'TII)2> 



for p' 



p-i- 



□ 



We will give a more general result of this type in Theorem \l.2\ 

5. The derivation operator 

In this section we define a derivation operator on our non- commutative proba- 
bility space and show that it satisfies similar properties as the derivation operator 
on Wiener space. 

We want to interpret the expression in the integration by parts formula in 



Proposition [4.7| as a directional or Frechet derivative. 
Definition 5.1. Let k ei)c® C^. We set 

DomDfc = i^B e B{Sj) -[Q{ki) — P{k2),B] defines a bounded operator on 

and for B e BomDk, we set DkB = f [Q(A;i) - P{k2),B]. 

Note that B E Dom Dk for some A; G [)c ® implies B* G Dom Dj: and 

DjB* = {DkB)*. 



Example 5.2. Let k G Ijc^C"^ and let G DomP(A;2)nDom(5(A;i)nDomP(A;2)n 
DomQ{ki) be a unit vector. We denote by the orthogonal projection onto the 
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one-dimensional subspacc spanned by ip. Evaluating the commutator [Q{ki) — 
P(/i;2),P^] on a vector e DomP(/c2) H DomQ{ki), we get 

[Q(h)-P(k2),W>^]<l> = {i;,<l>){Q(k,)-P{k2)){i^)-{i^,{Q(k,)-P(k2))(m 

= {^|;,c|>){Q{k,)-P{k,)){^|;) - {{Q{k,) - P{k2))^l^A)i^ 

We see that the range of [Q{ki) — P(A;2),P^] is two-dimensional, so it can be 
extended to a bounded operator on S^. Therefore e Dom£)fe, and we get 

(P>,Pv,)0= -P(A;2))W - {{Q{ki) - P{k2))i^A)i) 

for all G ij. 



Example 5.3. Let /i e ® E^^ A; G ()c ® C^. Then ^[Q{ki) - P{k2),U{hi, /i2)] 
defines a bounded operator on Sj, and we get 

DkU{hi, h2) = hi) + {k2, h2))U{hi, /12). 



Proposition 5.4. Let k e i)c<S>'C^- The operator is a closable operator from 
B{Sj) to B{Sj) with respect to the weak topology. 

Proof. Let (i?n)neN Q DomD^ C B{Sj) be any sequence such that Bn ^ and 
DkBn — > P for some /3 e B{S)) in the weak topology. To show that Dk is 
closable, we have to show that this implies P = 0. Let us evaluate (3 between two 
exponential vectors S{hi), S{h2), hi,h2 G i)c, then we get 

{£(hi),(3£(h2)) = \im{S(hi),DkBnS(h2)) 

n— >oo 

= lim '-{{Qih) - P(k2))S{hi),B^£{h2)) 

- lim U8{hi),B,,{Q{ki) - P{k2))S{h2)) 
= 0, 

and therefore /? = 0, as desired. □ 



Definition 5.5. We set 



5 = alg Oh{v) 



s.t. ^ ^ e DomOh for all Ki,K2>0 



dx'^^ dy 



«2 



the elements of S will play the role of the smooth functionals. Note that S 
is weakly dense in B{S)), i.e. S" — B{S^), since S contains the Weyl operators 
U{hi,h2) with /ii,/i2 e t)- 
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We define D : S ^ ^{^) ® f)c ® (where the tensor product is the algebraic 
tensor product over C) by setting DOh{<f) equal to 



( 



\ 



On 



dy 



® h2 



J 



and extending it as a derivation w.r.t. the B(io)-bimodule structure of B{Sj) 
f)c ® defined by 



O- 



01 (g) ki 

02 k2 



001 ® ki 

002 ® k2 



01 (g) ki 

02 ® k2 



■O 



OiO (g) ki 
O2O ® k2 



for O, Oi, O2 e B{S^) and A; e f)c C^. 



Example 5.6. For /i e f) (g M^, we get 

= i5o.(expi(.+,))=,([;(;;;;;;j|;:;) 

= iU{hi,h2)®h. 

Definition 5.7. We can define a -B(^)-valued inner product on B{^) (g) l^c ® 
by (•, •) : E{^) g) «) C2 X B{f)) ®\)c®C^ ^ B{S)) by 

We have 

= {A,Br 

0*{A,B) = (AO,S) 
(A,S)0 = {A, BO) 
{0*A,B) = (AOS) 

for all A,B e B{S)) ® f)c and all O e This turns B{^) ® i)c ^ 

into a pre-Hilbert module over -6(^5). It can be embedded in the Hilbert module 
Tl = B{S),S) (g [)c ® C^) by mapping O ®k E B{S^) (g [)c ® to the linear map 
ij^vi— >-Ow(gA;G^(gf)c<gC^. We will regard f)c g) as a subspace of 971 via 
the embedding l)c 9 /c >— >^ id^, (g /c G SDT. Note that we have O-k — k-O — 0®k 
and (A, k) = Z) for all k e\)c®'C^,0 e B{Sj), AeWl, where the conjugation 
in Wl is defined by UWk ^ O* ®k. 

Proposition 5.8. Let O E S and /c e f)c ® C^- Then O e DomD^ and 



DkO^{k,DO) = {DO,k). 
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Proof. For /i G (g) and G Bom Oh s.t. also ^ G Bom Oh, we get 



{k,DOh{v)) 



k2 



dx ' dy 



\ 



o, 



dip 
dy 



-[Q{h)-P{h),Oh{v)]=DkO, 



where we used Proposition [4.7| . The first equahty of the proposition now follows, 
since both O ^— DkO = ^[Q{ki) — P{k2),0] and O t-^ {k, DO) are derivations. 
The second equality follows immediately. □ 



The next result is the analogue of Equation (|2] 
Theorem 5.9. We have 

E{(k, DO)) = ^E{{P{k,) + Q{k2), O}) 

for all G f)c ® o.iT'd all O E S, where {■,■} denotes the anti- commutator 
{X,Y} = XY + YX. 

Proof. This formula is a consequence of the fact that Q{h)Q = h = iP{h)Q for 
all h G i)c, we get 

i 
2 



E{{k,DO)) 



ki) - P{k2))n, on) - {n, o(g(fci) - p{h 



-{{h + ih, OQ) - {Q, 0{ki + ik2))) 
]^[{{P(kr)+Q{k2))n,0n) + {n,0{P{kr) + Q{k2 

\mp{k{)+Q{k2),o}). 



□ 



There is also an analogue of 
Corollary 5.10. Let A; G t)c ® C^, and Oi, . . . , 0„ G 5, then 

-E ( <j p{kr) + Q{k2), n M = ^ E n ^^-) n 

) / \m=l j=l j=m+l 



m=l 



where the products are ordered such that the indices increase from the left to the 
right. 



Proof. This is obvious, since O t-^ {k, DO) is a derivation. 



□ 
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This formula for n = 3 can be used to show that D is a closable operator from 

B{s^) to mt. 

Corollary 5.11. The derivation operator D is a closable operator from B{S^) to 
the B{S))-Hilhert module 0?t = B{S^,S^ ®^ic® C^) w.r.t. the weak topologies. 

Proof. Wc have to show that for any sequence in S with A„ and 

DA,, a G 971, we_get a = 0. Let f,g e Set /i = /a = 

^1 = ^, and 52 = ^, then we have Uif^J^)^ = e~\\f\\'^8{f) and Uig^.g^)^ = 
e-M\l^8{g). Thus we get 

e^v{{\\f\? + \\g\?)/2){£{f)®h,a£{g)) 
= exp((||/|p + ||5|p)/2)(^(/),(7^,«)^(5)) 
= hm E(f/(-/i,-/2)(^,DA„)^7(5i,52)) 

n— +00 

= hm e(^{P(/ii)+Q(/i2),C/(-/i,-/2KC/(^i,^2)} 

-{h^DU{-f^,-f2))A^U{g,,g2)-U{-f,,-f2)A^{h,DU{g,,g2))) 

= lim ((-01, Ani^2) + (V'3, ^-04) - (-05, ^nV'e) " (-07, ^^nV's)) 
n—*oo 

= 

for all /i e l^c <8) C^, where 
V'l = |C/(/i,/2)(^(/^i) + Q(M)^^, V'2 = U{g,,g2)n, 

V^5 = (l)^f/(-/i,-/2))*f^, V-e = f/(^?i,^?2)l^, 

V-T = U{frj2)n, ^8 = D,f/((?i,(72)^^. 

But this implies a = 0, since {S{f) ® G f)C)^ & i)c ^ C^} is dense in 

Remark 5.12. This implies that D is also closable in stronger topologies, such as, 
e.g., the norm topology and the strong topology. 

We will denote the closure of D again by the same symbol. 

Proposition 5.13. Let O e DomD. Then O* G DomD and 



DO* = DO. 

In particular, since D is a derivation, this implies that DomI? is a *-subalgebra 
ofBiSj). 

Proof. It is not difficult to check this directly on the Weyl operators U{hi,h2), 
hei)(^R'^.We get U{hi, ^2)* = U{-hi, -h2) and 

D{U{hi,h2T) = DU{-hi,-h2) = -iU{-hi,-h2)®h 
= U{hi,h2y ®W) = DU{hi,h2). 
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By linearity and continuity it therefore extends to all of DomD. □ 

We will now show how D can be iterated. Let if be a complex Hilbert space, 

then wc can define a derivation operator D : S (E) H ^ ^{^) ^ i)c ® (E) H 
by setting D(0 h) = DO h ior O G S and h E H. Closing it, we get 
an unbounded derivation from the Hilbert module B{S^,S^ ® H) to ^M{H) — 
B{S) (g) if, ^ (8) {)c ® (8) if). This allows us to iterate D. It is easy to see that D 
maps S(^H to S^ijc'^C'^^H and so we have L>"(5®ii) C 5® ([)c ® C2)®''®ii. 
In particular, S C DomZ^" for all n G N, and we can define Sobolev-type norms 
1 1 • 1 1„ and semi-norms 1 1 • 1 on S by 

\\0\\l = \\0*0\\ + j2\\{D''0,D-0)\\, 

n 

\\o\\i,n = \m\'+j2m,{D''o,D-o)m for^eiD 

In this way we can define Sobolev-type topologies on DomD". 

We will now extend the definition of the "Prechet derivation" to the case 
where k is replaced by an element of 971. It becomes now important to distinguish 
between a right and a left "derivation operator" . Furthermore, it is no longer a 
derivation. 

Definition 5.14. Let m G 3Jl and O G DomD. Then we define the right gradient 
D uO and the left gradient O D u ol O with respect to u by 

'^uO - {%D0), 

o'd^ = {DO,u). 

We list several properties of the gradient. 
Proposition 5.15. 1. Let X G B{S^), O, Oi, O2 G DomL>, and ue^Jl. Then 

^«(0l02) = (^uOl^02 + ^uO,02, 

O^Dux = {0"d^)X, 
{Oi02)'d^ = Oi'd + Oi , 

2. For A; G t)c ® and O G DomZ^, we have 

DkO = Didf^i^kO = ODidf,m 

Proof. These properties can be deduced easily from the definition of the gradient 
and the properties of the derivation operator D and the inner product (, ). □ 
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Remark 5.16. We can also define a two-sided gradient D „ : DomL> x DomD 
B{S)) by ■■ (Oi,02) ^ 0,"^^02 = Oi l'^u02) + {pi"Du) O2. For k e 
® we have Oi^^id^j^fcOa = Dfc(0i02). 

6. The divergence operator 

The algebra B{SS) of bounded operators on the symmetric Fock space $^ and 
the Hilbert module OJl are not Hilbert spaces with respect to the expectation in 
the vacuum vector Q. Therefore we can not define the divergence operator or 
Skorohod integral 5 as the adjoint of the derivation D. It might be tempting to 
try to define 5X as an operator such that the condition 

(6.1) W.{{5X)B) = E {^xB^ 

is satisfied for all B G Dom D x, even though it is not sufficient to characterize 
5X. But the following proposition shows that this is not possible. 

Proposition 6.1. Let /c e {)c ® with ki + ik2 7^ 0. There exists no (possibly 
unbounded) operator M whose domain contains the vacuum vector such that 

E{MB) =E{DkB) 

holds for all B G DomDk. 

Proof. We assume that such an operator M exists and show that this leads to a 
contradiction. 

Let B e B{S^) be the operator defined by Sj 3 ip ^ {ki + ik2, i^)^, it is easy to 
see that B e Dom and that D^B is given by 

{DuB)ij ='^{ki + tk2, ip){ki + ik2) - '-{{Q{ki) - P(k2)){k, + zfcs), ^>fi, for ^ e i^- 
Therefore, if M existed, we would have 

= (O, MBn) = E{MB) = E{DkB) 

= {n,{DkB)n)^~{ki + ik2,ki + ik2), 
which is clearly impossible. □ 
We introduce the analogue of smooth elementary [)- valued random variables. 

If we define AtuB for u = E"=i Fj (g) h^^ e St, and A, Be B{Sj) by 

n n 
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then it follows from Corollary p.lO| that this satisfies 
(6.2) E (^AtuB^ = E (^A^uB^ . 

But this can only be written as a product AXE with some operator X, if A and 

B commute with P (h^^^ + Q (li:^^^ , . . . , P (hf^^ + Q (hf^^ . We see that a 

condition of the form (|6.1|) or (|6.2|) is too strong, if we require it to be satisfied 
for all A,Be DomD. We have to impose some commutativity on A and B to 
weaken the condition, in order to be able to satisfy it. We will now give a first 
definition of a divergence operator that satisfies a weaker version of ( |6.2|) , see 



Proposition |6.3| below. In Remark |6.6| we will extend this definition to a bigger 
domain. 

Definition 6.2. We set 

n -. n n 

j=i j=i j=i 

defines a bounded operator on S) 
and define the divergence operator 6 : St^^s B{S)) by 

^ n n 

for u = Fj ® h^'^ e S^^s- 

It is easy to check that 

m) = {siu)y 

holds for all u G Stj^s. 

Proposition 6.3. Let u = ^J^i Fj ® h^^"* G Si^^s 

A, 5 G DomD n |P (/iS'^) + Q {h^^^^ , . . . , P (/li")) + Q (^"^ 

i.e., A and B are in the commutant o/|p (h^^'^^ + Q {h^^'^^ , . . . , P (h^i^^ + Q (h 
then we have 

E(v45(m)P) = E [a'TSuB 

Remark 6.4. Note that 5 : S^^s 'S(^) is the only linear map with this property, 
since for one single element /i G f)c ® C^, the sets 



A*n 

BVL 



A G DomP) n {P {hi) + Q (/is)}' 
B G DomP n {P {hi) + Q (/is)}' 
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are still total in S^. 

Proof. From Corollary |5.10| we get 

E ( A^„B 



E{A{u,DB) + {DA,ti)B) 

(n n 
i=i i=i / 



But since A and B commute with P {h^^^ + Q {h^^'^^ ,P (h^^^^ + Q (h 
we can pull them out of the anti-commutator, and we get 



2 h 



E AD„B 



n n ^ 

- 5^ A {P [h?) + Q (/.(^■)) ,F,]b-Y.A (D.u^F,) B 

j=l j=l y 



E{AS{u)B). 



□ 



We will now give an explicit formula for the matrix elements between two 
exponential vectors of the divergence of a smooth elementary element u G Stj^s, 
this is the analogue of the first fundamental lemma in the Hudson-Parthasarathy 
calculus, see, e.g., ||Par92| , Proposition 25.1]. 



Theorem 6.5. Let u G S^^s- Then we have the following formula 

iki — ik2 



{S{k^),6{u)£ik2)) = {£{k,) 



ki + k2 



, u£{k2 



for the evaluation of the divergence S{u) of u between two exponential vectors 
8{ki), £{k2), for ki, k2 G i)c- 

Remark 6.6. This suggests to extend the definition of 6 in the following way: set 



Dom5 = |m G OJl 3M G ?B(.^),VA;i,A;2 G f)c : 
(6.3) {£{k,),M£{k2)) = (^S{k,) ® ( ) ,^^(^2)^ } 

and define 6{u) for u G Dom6 to be the unique operator M that satisfies the 



condition in Equation (|0|). 
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Proof. Let u = Yl^=i Fj0h^^\ Recalling the definition of we get the following 
alternative expression for 6{u), 

n 

(6.4) = J2 (^2'^ - ^^S'O ^-^^ + (^^'^ ~ '^^'0 ) ■ 

Evaluating this between two exponential vectors, we obtain 

n 



n 



□ 



Corollary 6.7. The divergence operator 6 is closable in the weak topology. 

Proof. Let {un)n£N be a sequence such that — > and — > € in the 
weak topology. Then we get 



{£{k,),p£ik2)) = \im{£ik,),SM£ik2)) 

iki — ik2 



n— »oo 





for all ki, k2 G fjc, and thus /3 = 0. □ 
We have the following analogues of Equations ( p. 3D and ( |2.5| ). 
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Proposition 6.8. Let u,v E Sff^s, F e S, h E l)c'Si C^, then we have 



(6.5) Dho5{u) = {h,u) + 5oDh{u) 

(6.6) 5iFu) = Fdiu)-F'Du + -J2[p{h^i'^)+Q{hi'^),F'\Fj 



(6.7) S{uF) = S{u)F-^uF + -J2[F,p(^h'f^)+Q(^hi''')]F, 
Proof. l.Letu^ ^Li Fj h^^\ We set 



Y 



-{P{h^^^]+Q{h, 



Q{hY']-P{h 



'2 ^ 



-{Q{h)-P{h)), 



then we have 6{u) = YTj=i {{^j - ^3)^3 + ^A^j + and therefore 



Dh{5{u)) = {y{X, - Y,)F, + YF,{X, + F,) - (X, - Y,)F,Y - F,(X, + Y,)Y). 
3=1 

On the other hand we have Dh{u) = ELi(^-P'j - FjY) ® h^^\ and 



s{D,{u)) = J2 {i^3 - y3)yF3 - i^3 - y3)F3y + yF3{x, + r,) - f,y{x, + r,)). 

Taking the difference of these two expressions, we get 



Dn{5{u)) - S{D,{u)) = ^ {[Y,X, - Y,]F, + + F,]) 
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2. A straightforward computation gives 



SiFu) = X^((X,-r,)FF,+FF,(X,+F,)) 



n n 
n n 
n 71 



where we used that [Xj, F] = i ^ , DF^ defines a bounded opera- 

tor, since F & S DoniD. Equation (|6.7|) can be shown similarly. 



□ 

If we impose additional commutativity conditions, which are always satisfied 
in the commutative case, then we get simpler formulas that are more similar to 
the classical ones. 

Corollary 6.9. If u = ^"^^ Fj ® h^^^ E and 

FeBomDn{p (/iS^)) + Q (4^^) , . . . , P (/iS")) + Q (4"^) }' 
then we have 

6{Fu) = F6{u)-F'd^, 
6{uF) = 5{u)F-l5uF. 

7. Examples and applications 

7.1. Relation to the commutative case. In this section we will show that 
the non-commutative calculus studied here contains the commutative calculus, 
at least if we restrict ourselves to bounded functionals. 

It is well-known that the symmetric Fock space r([)c) is isomorphic to the 
complexification L^(f2;C) of the Wiener space L^(fi) over f), cf. ||Bia93| , |Jan97 , 

|Mey95|] . Such an isomorphism / : L^(fi; C) ^ r(f)c) can be defined by extending 
the map 

/ : ^ J ^^^w(h)^^ ^ ^iQih)^ ^ e-ll^ll'/2£:(z/i), for he^. 
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Using this isomorphism, a bounded functional F G C) becomes a bounded 

operator M(F) on r([)c), acting simply by multiplication, 

M{F)^ = l{FI-\iP)), for G r(f)c). 

In particular, we get M [e^^i^)) = f/(0, h) for /i G 1). 

We can show that the derivation of a bounded differentiable functional coin- 
cides with its derivation as a bounded operator. 

Proposition 7.1. Let k e and F e C) nDomD^ s.t. DkF G L°°(^];C). 

Then we have M{F) G DomD^Q, where = ^ ^ ^ ; one? 

MiD,F) = D,,{M{F)). 

Proof. It is sufficient to check this for functionals of the form F = e^^^'^\ /i G f). 
We get 

M{Dke'^^^^) = M {i{k,h)e'^^''^) 

= t{k,h)U{0,h)=z(( l))uiO,h) 



^ ^k J ' \ h ^ ^ 

□ 

This implies that we also have an analogous result for the divergence. 

7.2. Sufficient conditions for the existence of smooth densities. In this 
section we will use the operator D to give sufficient conditions for the existence 
and smoothness of densities for operators on Sj. The first result is a generalisation 
of Proposition ^4.8| to arbitrary states. 

Theorem 7.2. Let k G N, /i G f) Cg> with {hi, h2) ^ 0, and suppose that $ is 
of the form 

^{X) = tr(pX) for all X G 
for some density matrix p. If there exist k,i & f)c ^ such that 

andpe n«^+^^<^DomD^^Df , and 

tT{\D'^Wpp\) < oo for all Ki + K2 < k, 

then Wh^^ G f]2<p<oo -^^''^(^^)' ^-^^ Wigner density Wh,<i> lies in the Sobolev 
spaces of order n. 
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A 



and set 
then we have 

[X2,0h{ip)] 



Xi\_^ ,-1 f Qih) - P{k2) 

"2^ Q{ii)-p{i2) y 



{h2,e2)DkOh{v) - {h2M)DiOh{v) _ 
del A 



detA 



for all Schwartz functions </?. Therefore 



/ 



QK1+K2 



= tr(p ]Xi, .. . [Xi , [X2, .. . [X2 , Ohif ) 

Ki times K2 times 
= \ti{[X2,...[X2,[Xi,...[Xi,p]]]]Oh{ip))\ 

< Cp,Kl,K2\\Oh{v)\\ < Cp,Kl,K2Ch,p\\v\\p, 

for all p e [1,2], since p e n«,+«,<« DomL'^iL'f and tT{\D'^' D^' p\) < 00 for all 
Ki + K2 < K, and thus 



C, 



P,K1,K2 



tr\[X2,...[X2,[Xi,...[Xi,p]]]]\<oc. 



But this implies that the density of dW^,* is contained in the Sobolev spaces 
HP^%R^) for all 2 < p < 00. ' □ 

Example 7.3. Let < Ai < A2 < ■ ■ ■ be an increasing sequence of positive 
numbers and {ej\j e N} a complete orthonormal system for l)c- Let : f)c — i)c 
be the contraction semigroup defined by 



for j en,t>o, 



with generator A = XljeN'^i^i- sequence increases fast enough to ensure 

that e~*'^J < 00, i.e. if trT^ < 00 for t > 0, then the second quantization 

Pt — T{Tt) : — > is a trace class operator with trace 

= trpt = ^ (en,Pten) 

neNf 
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where we use NJ^ to denote the finite sequences of non-negative integers and 
{en|n G Nj°} is the complete orthonormal system of Sj consisting of the vectors 

e„ = ef"i © ■ ■ ■ ef for n = (rii, . . . , n,.) G Nf, 

i.e. the symmetrization of the tensor ei ® • ■ ■ ei ® ■ ■ ■ (S> ©> ■ ■ ■ where each 
vector Cj appears Uj times. We get Zt = XlneN- UkLi 6""'=*^'= = [Ifcli i-e-^'^k 
the trace of pt. We shall be interested in the state defined by 

<I>(X) = ^tr(piX), for X G B{Sj). 



We get 

XI |(en, |pi/2«^(ej)ren)| = ^ \\pti2a\tj) 



|2 



00 ^ 



n=0 /cT^j 

and therefore that pta^{ej) defines a bounded operator with finite trace for all 
j, £ G N and t > 0. Similarly we get 

tr I a^(ej)pt I < 00, tr |pt (a^(ej))^| < 00, etc. 

and 

ti \P^' {ej^)Q^^ (ej^) pt\ < 00, tr |P^i(ej-JpiQ^2(ejJ| < 00, etc. 

for alH > and ji, J2, ^1, ^2 e N. For a given G P) ® with {hi, ^ (and 
thus in particular hi and /i2 7^ 0), we can always find indices ji and j2 such 
that {hi,ej^) 7^ and {h2.,ej^) 7^ 0. Therefore it is not difficult to check that for 



all K G N all the conditions of Theorem |7^ are satisfied with k = { | and 







^32 



We see that the Wigner density Wh,'^ of any pair {P{hi),Q{h2)) 



with {hi,h2) ^ in the state $(■) = tT{pf)/Zt is in fl.eN n2<p<oo 

in particular, its derivatives of all orders exist, and are bounded and square- 

integrable. 

We will now show that this approach can also be applied to get sufficient con- 
ditions for the regularity of a single bounded self-adjoint operator, for simplicity 
we consider only vector states. 
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Given a bounded self-adjoint operator X, we call a measure /ix,* on the real 
line its distribution in the state if all moments agree, 



= [ x"d/ix,4., for all n G N. 
Jr 



Such a measure always exists, it is unique and supported on the interval 

[-II^IUI^I 



Proposition 7.4. Let u & S) be a unit vector and let $(■) = {uj,-uj) be the 
corresponding vector state. The distribution nx,<s> of an operator X G B{S}) in 
the state $ has a bounded density, if there exists a k E i)c ® such that uj G 
Dom(Q(A;i)-P(A;2)) nDom(g(A?i) -P(^2)), X G DomD^, X-DtX = DtX-X, 
D^X invertible and {D^X)^^ G DomZ^^. 

Proof Since X ■ D^X = D^X ■ X, we have DkPiX) = {DkX)p'{X) for all poly- 
nomials p. We therefore get 

{{DkX)-'p{X)) =p{X)Dk {{DkX)-') +p'{X). 

The hypotheses of the proposition assure 

\{uj,Dk[{DkX) p{X))lj)\ < ^ \\{DkX) ||||p(X)|| 

< Ci sup |p(2^)L 

^6[-||X||,||X||] 

\{u,piX)D,{iDkX)-')co)\ < \\D {iD,X)-')\\ MX)\\ 

< C2 sup \pix)\, 

x&[-\\X\\,\\X\\] 

and therefore allow us to get the estimate 
^11 



p'{x)dfix,<s>{^) 



\{lu,p'{X)lo)\ 



00, -p{X)Dk {{DkX)-') ju 

< (Ci + C2) sup 1^(2;) I 

^6h||X||,||X||] 

for all polynomials p. But this implies that admits a bounded density. □ 

Let n G N, n > 2. We get that the density is even n — 1 times differentiable, 
if in addition to the conditions of Proposition O we also have X G DomZ?^, 
{DkXy^ G DomD^, and 

uje f] Dom{Q{h) - P{k2)y fl Dom(g(fci)-P(fc2))". 

l<K<n l<K<n 
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The proof is similar to the proof of Proposition [7.4| , we now use the formula 

where Aq, . . . , An-i are some bounded operators, to get the necessary estimate 

< C sup 



\\x\\ 



\X\\MX\ 



\p{x)\ 



by induction over n. 



7.3. The white noise case. Let now f) = L^{T,B, fx), where {T,B,fi) is a mea- 
sure space such that B is countably generated. In this case we can apply the 
divergence operator to processes indexed by T, i.e. i3(^)-valued measurable func- 
tions on T, since they can be interpreted as elements of the Hilbert module, 
if they are square-integrable. Let L'^(T,B{S))) denote all -valued measur- 
able functions t ^ Xf on T with Jj, ||X(|pdt < oo. Then the definition of the 
divergence operator becomes 

Dom5 = 



X = {X\X^) G L^T,B{^))(BL^T,B{^)) 
3M G B{Sj),yh,k2 G f)c : {S{h),MS{k2)) 

{t{h-h){£{h),X\t)£{h)) + {h + k2){£{h),X\t)£{h 



and S{X) is equal to the unique operator satisfying the above condition. We will 
also use the notation 



5{X) 



X\t)dP{t) + [ X\t)dQ{t), 
Jt 



and call S{X) the Hitsuda-Skorohod integral of X. 

Belavkin [|Bel91a| , |Bel91b|] and Lindsay [[Lin93|] have defined non-causal quan- 
tum stochastic integrals with respect to the creation, annihilation, and conser- 
vation processes on the boson Fock space over L^(M+) using the classical deriva- 
tion and divergence operators. It turns out that our Hitsuda-Skorohod integral 
coincides with their creation and annihilation integrals, up to a coordinate trans- 
formation. This immediately implies that for adapted, integrable processes our 
integral coincides with the quantum stochastic creation and annihilation integrals 
defined by Hudson and Parthasarathy, cf. [ 1IP84 , |Par92 ]. 

Let us briefly recall, how they define the creation and annihilation integral. 



cf. [[Lin93|| . They use the derivation operator D and the divergence operator 6 
from the classical calculus on the Wiener space L^(f2), defined on the Fock space 
r(L2(M+; C)) over L'^(R+; C) = L2(M+)c via^the isomorphism between L'^{n) and 
r(L^(]R+; C)). On the exponential vectors D acts as 

D£{k) = £{k) ® k, for A; G L2(M+,C), 
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and 5 is its adjoint. Note that due to the isomorphism between r(L^(R+; C)) 
(R+ ; C) and (^R+ ; T (L^ (R+ ; C) ) ) , the elements of T (L^ (R+ ; C) ) (g) (R+ ; C) 

can be interpreted as function on R+. In particular, for the exponential vectors 
we get (^DS{k)^^ — k{t)£{k) almost surely. The action of the annihilation integral 
J FfdAt on some vector -0 G r(L^(R+; C)) is then defined as the Bochner integral 

/•BL /. 

and that of the creation integral as 

BL 

FtdA^il) = ~5{F.^P). 

■+ 

These definitions satisfy the adjoint relations 

(/•BL \ * /.BL / /.BL \ * ^BL 

/ FtdAr] D / F^dAl and / FtdA;] D / F^dA^. 

Proposition 7.5. Let {T,B, ji) = (R+, i3(M+), da;), i.e. the positive half-line with 
the Lebesgue measure, and let X — {X^,X^) e Dom5. Then we have 

I- fBL /-BL 

X\t)dP{t)+ / X\t)dQ{t)^ / {X^-iX^)dA;+ / {X^ + iX^)dAt. 

Proof. To prove this, we show that the Belavkin-Lindsay integrals satisfy the same 
formula for the matrix elements between exponential vectors. Let {Ft)te'M.+ £ 
L^(R+,B(i3)) be such that its creation integral in the sense of Belavkin and 
Lindsay is defined with a domain containing the exponential vectors. Then we 
get 

/■BL 

{£ih), FtdAlSih)) = {£{h),S{F£{k2))) 

= {{DS{k,)),FS{k2)) 



k,it){£ik,),Ft£{k2))dt. 



For the annihilation integral we deduce the formula 

r.BL /-BL 



{£{k,), F,dA£{k,)) = (/ F:dA^£{k^),£{k2)) 



= / k2{t){£{k2),F*£{ki))dt 

Jr^ 

k2{t){£{k,),Ft£{k2))dt. 
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□ 

The integrals defined by Belavkin and Lindsay are an extension of those defined 
by Hudson and Parthasarathy, therefore the following is now obvious. 

Corollary 7.6. For adapted processes X G Dom6, the Hitsuda-Skorohod integral 

6iX)= [ X\t)dPit)+ [ X\t)dQit) 
Jt Jt 

coincides with the Hudson- Parthasarathy quantum stochastic integral defined in 



7.4. Iterated integrals. We give a short, informal discussion of iterated inte- 
grals of deterministic functions, showing a close relation between these iterated 
integrals and the so-called Wick product or normal-ordered product. Doing so, 
we will encounter unbounded operators, so that strictly speaking we have not 
defined 5 for them. But everything could be made rigorous by choosing an ap- 
propriate common invariant domain for these operators, e.g., vectors with a finite 
chaos decomposition. 

In order to be able to iterate the divergence operator, we define 5 on B{SS) ® 
f)c ® ® if, where H is some Hilbert space, as 5 ® idj^. 

Using Equation (|6.4|), on can show by induction 



2 / V "2 / /C{l,...,n} je/ ie{l,...,n}\/ 




in) 



e f)c ® This is just the Wick 
.,(P(/^S"^ + g(/.("))),i.e. 

= {Pih['^ + Q{hP)) o . . . o (P(/^S") + Q{hP)) , 

where the Wick product o is defined on the algebra generated by {P{k), Q{k)\k G 
f)c} by 

P{h)oX = X oP{h) = -ia+{h)X + iXa(h) 
Q{h)oX = X oQ{h) = a+{h)X + Xa(h) 

for X G alg{P{k),Q{k)\k G f)c} and /i G f)c in terms of the momentum and 
position operators, or, equivalently, by 

a-^{h)oX = X oa+{h) = a+{h)X 
a{h)oX = Xoa{h) = Xa{h) 

in terms of creation and annihilation. 
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8. Conclusion 

We have defined a derivation operator D and a divergence operator 5 ojiB{SS) 
and S)®\)c® C^), resp., and shown that they have similar properties as the 
derivation operator and the divergence operator in classical Malliavin calculus. 
As far as we know, this is the first time that D and 5 are considered as operators 
defined on a non-commutative operator algebra, except for the free case |[BS98 



where the operator algebra is isomorphic to the full Fock space. To obtain close 
analogues of the classical relations involving the divergence operator, we needed 
to impose additional commutativity conditions, but Proposition |6.1| shows that 
this can not be avoided. Also, its domain is rather small, because we require 5{u) 
to be a bounded operator and so, e.g., the "deterministic" elements /?. G [)c ® 
are not integrable unless h = 0. One of the main goals of our approach is the study 
of Wigner densities as joint densities of non-commutating random variables. We 
showed that the derivation operator can be used to obtain sufficient conditions 
for its regularity, see, e.g.. Theorem [7.2| . It seems likely that these results can be 
generalized by weakening or modifying the hypotheses. It would be interesting 
to apply these methods to quantum stochastic differential equations. 
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